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The 2D system of electron confined to the lowest Landau level is described using a representation of 
the density matrix depending both on electron and hole coordinates. Condensation of the electron 
system into a fractional quantum Hall state is associated to clustering of particle and hole coordi- 
nates. The correlation of particle and hole coordinates is studied and groundstate wave functions 
are derived for v = p/(2p ± 1). These wave functions prove to be accurate for the studied cases 
v — 2/3 and 3/5 and are identical to pair wave functions [R. Morf et ah, Phys. Rev. B34, 3037 
(1986)] for v = 2/3 and 2/5. 
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The theory of Laughlinlil provides a good microscopic 
description of the v = 1 /to incompressible quantum Hall 
states, with to odd. Based on this description, it has been 
proposed to characterize fractional quantum Hall states 
by the occurence of binding of flux quanta to electrons!! 
An explicit formulation for such a .-binding is provided 
by the composite Fermion approacbtra, which makes use 
of a singular gauge transformation to map the system of 
electrons in presence of magnetic field onto a system of 
composite Fermions, which are objects composed of an 
electron bound to an even number of flux quanta. 

In this note, an alternative formulation is presented 
for treating the binding of flux quanta to electrons. This 
approach makes use of a particular representation of the 
finite temperature density matrix for the electron system 
confined to the lowest Landau level which consists in ma- 
trix elements of the density matrix operator evaluated 
between states parametrized by electron coordinates and 
states parametrized by the coordinates of electron vacan- 
cies, i.e. holes. The latter coordinates can be regarded 
as positions of flux quanta. 

This formulation, in which electron-hole symmetry is 
fully explicit, leads to a simple and natural picture for 
the incompressible quantum Hall state at v = p/q: At 
some finite temperature, the density matrix is large only 
for configurations in which electron and hole coordinates 
can be grouped into clusters containing each p electron 
coordinates and q — p hole coordinates. This electron- 
hole binding results dynamically from the repulsion be- 
tween electrons. The creation of a neutral excitation cor- 
responds to the breaking of such a cluster into two smaller 
clusters, each one leading to a local fractional charge de- 
fect. The fact that this breaking requires a finite amount 
of potential energy leads to a gap in the energy spectrum 
and to incompressibility. The off-diagonal long range or- 
der characterizing incompressible groundstatesEl results 
from the fact that these clusters are Bosons and Bose- 
condense. 

The advantage of working with this finite temperature 
description is that the correlation between hole coordi- 
nates, which play the role of positions of flux quanta, 



and electron coordinates can be written down explicitely. 
However, the price to pay is the work necessary to de- 
rive from this density matrix a microscopic description 
of the system a zero temperature. We show how to 
obtain groundstate wave functions either as a function 
of electron coordinates or as a function of hole coordi- 
nates in the case of the fractional quantum Hall states 
at v = p/ (2p +1)- The case is made that the wavefunc- 
tions derived from this density matrix are generalizations 
of the pair wavefunctions of Ref.tl, which are known to 
be accurate. In the latter wave functions, the pairing of 
electrons can be regarded as resulting directly from the 
presence in the density matrix of clusters containing each 
p = 2 electrons and q — p holes. 

This description of fractional quantum Hall states 
based on a density matrix is in fact analogous .-to the 
one provided by the composite Fermion approachou. In- 
deed, a cluster made up of one electron and one hole 
corresponds to a composite Fermion and a cluster made 
up of p electrons and p + 1 holes, which enters into the 
description of the v — p/(2p+ 1) state, corresponds to 
p composite Fermions, each one belonging to a different 
Landau level, all bound to a flux quantum of the field 
experienced by composite Fermions. 

Let us first study the density matrix of the 2D elec- 
tron gas confined to the lowest Landau level at a filling 
factor v = p/q characterized by the presence of an in- 
compressible groundstate. More specifically, we consider 
the quantity 

p l3 (wi...w hN z 1 ...z pN ) = (w 1 ...w hN \e~ f3V \zi...z pN ) (1) 

where V is the electron-electron interaction projected 
onto the lowest Landau level and where \z\...z v n) and 
\w\...WhN) are states respectively given in electron rep- 
resentation and hole representation, defined by 

\z\...z pN ) = c)(zi)..J(z pN ) |0) , (2) 
\wi...w hN ) = c(w 1 )...c(w hN ) |1) , (3) 

where |0) and |1) respectively denote the empty- and full- 
Landau level states, c^(z) is the operator creating an elec- 
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tron in a coherent stateQ centered at position z, his given 
by h = q — p and where pN and hN are respectively the 
number of electrons and holes in the lowest Landau level. 

For (3 = 0, the density matrix of Eq. (|l|) is simply 
given by 
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in units of the magnetic length, where Zi and wj are re- 
spectively the complex coordinates of electron i and hole 
j. Ignoring for a moment the distinction between vari- 
ables Zi and Wj , we see that po has a functional form iden- 
tical to that of the wave function describing a full Lan- 
dau level. Thus configurations {w\...WhN z\...z p n) char- 
acterized by non-negligeable values of po have for feature 
( denoted by (A) ) that their coordinates are homoge- 
neously distributed over the whole sample. Let us now 
study the effect of the electron-electron repulsion present 
in Eq. ([!]). In the case of a repulsiorj|-|Characterized by 
vanishing pseudo-potential coefficients^ for even values 
of the electron-electron relative angular momentum, the 
Hamiltonian matrix is given, up to an additive constant, 
by 



d 



pN hN 



-QpP\p=o(wi...w hN z 1 ...z pN ) = -xmH^ 



x / d 2 ^ 



i=l 3 = 1 L 

{wi...w hN z 1 ...z pN ) 1 (5) 



where the sum over L is carried over odd positive integers, 
Vl denote the pseudo-potential coefficients, ip^ ( z i z z) is 
the wave function describing two particles with relative 
angular momentum L and with center of mass in a co- 
herent state centered at position £ and where V is the 
projector associated to this wave function, defined by 



{wi...w hN z 1 ...z pN ) = ip]p(ziWj) x 



d 2 z' i d 2 w' j ip < ^' (z' i w' j )p (wi..w' j ..w hN z 1 ..z' i ..z pN ). (6) 



As can be seen in Eq. (g), the effect of the interaction 
at finite values of (3 is to increase pp for configurations 
which have for feature ( denoted by (B) ) that their co- 



ordinates 



are close to their coordinates 



Thus it 



is likely that the quantity pp for small but finite (3 will 
be largest for configurations {w\...WhN Z\...z p n) in the 2D 
plane which can be divided into compact clusters of equal 
sizes containing each p positions z% and h positions Wj, 
which possess both features (A) and (B). 

Let us now construct a form for the density matrix 
satisfying the latter requierement for the special case of 
the filling fractions v = p/(2p + I). The question of 



the temperature range in which this form is valid will be 
discussed latter. We propose to use the function given by 
n(w 1 ...w hN zi...z pN ) = A w A z k{wi...w hN zi...z pN ) and by 



f N 



X (5n) (Zp(n-l)+l---Zpn.H>/,(n-l)+l-" t»»n) 



Po, (7) 



where A w and A z are antisymmetrization operators 
acting respectively on sets of variables {wi. ..w^n) and 
{z\...z p n) and where x^"" 1 i s a function describing the 
correlation of p electrons and h = p + 1 holes in a cluster- 
centered at position £„, given by 



X ( - i) (z 1 ..z p w 1 ..w h ) = Y[ On ~ z i 2 ) II ~ w hf 
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The projectors and the term po hi Eq. (Q) respectively 
ensure that configurations leading to a large value of k 
have features (B) and (A). The integrals implicit in the 
projectors of Eq. ([?]) can be evaluated using the identity 
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PZ{z 1 ..z n )Qi{z 1 ..z n ), (9) 



where / is a polynomial and where Pk and Qi are ho- 
mogeneous polynomials of total degree k and I with 
k < I and, are invariant under global translations of their 
variablesfl This leads, up to a multiplicative constant, 
to 

r N 

k = Yl [ d2 ^X (in) ( z p(n-l) + l--- Z pnWh(n-l) + l---Whn) 
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The function x is choosen so as to satisfy the following 
criteria: (a) is antisymmetric under interchange of 
two electron coordinates or under interchange of two hole 
coordinates, (b) The total degree of the polynomial part 
of x 1 - ^ equals q(q — l)/2. Indeed, in light of Eq. (||), k 
would vanish for a smaller degree. On the other hand, a 
larger degree would leadu to a reduction of the power of 
factors (£ ni — £„ 2 ) in Eq. ( |l0| ) and configurations with a 
large value of k would not satisfy condition (A) anymore, 
(c) Factors of type (zi — Wj) are absent so that x is large 
for configurations satisfying condition (B). 

In order to support our claim that the function k pro- 
vides a realistic description of the density matrix for some 
temperature, we now make the case that the groundstate 
wave functions derived from this form for the density 
matrix are generalizations of the pair wave functions of 
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Ref.B, which are known to be accurate. Let us denote by 
$(zi...z p Ar) and Q(wi...WhN) the wave functions describ- 
ing the groundstate at v — pj (2p + 1) in the electron- and 
hole-representations of Eqs. (|2^ and (||), respectively. If 
k provides a good description of the density matrix, then 
$ and should maximize the quantity 

me)= f d 2pN Zd 2hN Wn*(W, Z)<fr(Z)Q(W) 
sj $ d 2 r N Z\$(Z)\ 2 j <P hN W \Q{W)f 

where Z = (zi...z p n) and W = (wi...WhN)- 

We now approximate -/?(<&, O) by substituting k for k 
in Eq. (|ll|). Ignoring for the moment the constraint of 
antisymmetrization, we consider wave functions of the 
type 

P N ^ N 

$(zi...Z pN ) = Y[e z — Y[ aijz (n _ 1)p+1 ...z np ) 

i—l n—1 

N 

X IT ^ l b( Z (n 1 -l)p+l--- z n 1 pi Z(„ 2 _ 1 ) p+1 ...Z n2 p) (12) 

n 1 <n 2 

where ai a and bi b are translationally invariant and ho- 
mogeneous polynomials 
of total degree l a and lb to be determined. The poly- 
nomial a; a (z( n _ 1 )p +1 ...z np ) describes correlations within 
the group n of variables zi n _is p+ i...z n p- The polyno- 
mial &z b (z(„ 1 _ 1 ) p+1 ...z„ 1 p,Z(„ 2 _ 1 ) p+1 ...2:„ 2 p) is built using 
products of l b factors of type (z( ni _ 1)p+n - z (n2 _ 1)p+i3 ), 
with ii,ia = 1— Pj and describes correlations between 
groups n i and n 2 . We write similarly 

Q(w!...W hN ) = *~ Y\_ C lai w (n-l)h+l--W n h) 

j=l n=l 

iV 

x ]~| rf; £i (w(„ 1 _i)/ l+ i...'!«„ 1 h,u)( n2 _i); 1+ i...w n2 / l ). (13) 

ni<«2 

We now introduce Eqs. ©, @, and (||) into Eq. 
( pd| ) and integrate over coordinates Zj and Wj before £ n 
in the numerator. In order to maximize the numerator, 
we demand that the result of the integration over elec- 
tron and hole coordinates have the same dependence in 
variables £ n as the kernel of k given by the second li 
of Eq. ( |Io| ) . Using Eq. (g) , it can be seen by inspectio: 
that this is possible only when l a and l c are equal to the 
degrees of polynomial parts of Eq. (||) depending respec- 
tively on coordinates z, and Wj, i.e. l a = p{p — l)/2 and 
l c = 3h(h — l)/2, and when l b + l d — q 2 . The latter con- 
dition, together with the constraint that each coordinate 
in $ or in O have qN — 1 zeros, leads to 4 = pq and 
Id = hq. In light of Eqs. (||) and (||), it is judicious to 
choose 



ai a (zi...z p ) = Y\_ ( 2 



(14) 



C lc {Wl...W h ) = f[ {Wj 



J 3 
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in order to obtain a large result from integration over 
particle and hole coordinates in the numerator of Eq. 
( |il"| ) . Among polynomials b pq and dh q leading to the same 
value for the numerator, those yielding to the most ho- 
mogeneous distribution of factors (z^ — z,- 2 ) over pairs 
1 < i\ < %2 < pN and of factors (wj 1 — Wj 2 ) over pairs 
1 < ji < j2 < hN lead to the smallest denominator. We 
therefore choose 

p p 

b pq {z 1 ...z 2p ) = n n ~ z p+i*f 

il = l 22=1 
P 

P P i=l 
/i ft 

d hq (wi...w 2h ) = n n ~ u, fc+ja) 2 

31=1 j 2 = l 
h 



and 



where P p and respectively denote sums over permu- 
tations of p and h objects. 

The wave functions $ = A z <t and = A w O de- 
fined using Eqs. ©, ©, (Q), ©, © agd @) 
are identical to the pair wave functions of Reffl when 
p = 2 and h = 2, respectively. $ and provide a good 
microscopic description of groundstates at filling factors 
v = p/(2p+ 1) and v = h/(2h — 1), respectively. Indeed, 
they lead to numbers of flux quanta which are in agre 
ment with the prescription of the hierarchical scheme! 
and have a large overlap with the exact groundstates 
of the corresponding systems of particles on a sphere in 
presence of Coulomb repulsion: In the case v — 2/3 and 
for 8 and 10 particles respectively, the overlaps are 0.954 
and 0.930 whereas the numbers of (L = 0) states are 2 
and 6. In the case v = 3/5 and for 6 and 9 particles 
respectively, the overlaps are 0.988 and 0.970 and the 
numbers of (L = 0) states are 3 and 8. 

We now study the constraints imposed by statistics on 
allowed values of parameters p, h and q = p + h using a 
path integral representation of the partition function 
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The exponential in Eq. (|T^) is brocken into M 
pieces exp(— (3V/M) with M even, and projectors 
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are inserted alternatively at imaginary times 
m0/M for m even and m odd, respectively. We con- 
sider configurations which can be divided at all imagi- 
nary times in clusters containing p coordinates zf"^ 1 and 
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h coordinates Wj , which significantly contribute to the 
partition function. For simplicity we further demand 

that z^ and Wj have a smooth dependence on m and 
belong to the same cluster throughout imagimary time. 
The phase corresponding to a given path is a product of 
phase contributions associated to each time step, which 
can in turn be evaluated by means of Eq. (0) . Ignoring 
factors (zi — Wj) for Zi and Wj belonging to the same 
cluster, this leads to (—l) ph ^ p \ where [P] is the parity 
of the permutation taking place among clusters between 
imaginary time r = and t — (3. Thus such paths con- 
tribute constructively to the partition function only if ph 
is even, that is if (a) p is even and h is odd or (b) p is 
odd and h is even or (c) both p and h are even. Note 
that cases (a) and (b) imply a filling factor v = p/q with 
odd denominator. 

The fact that all incompressible fractional quantum 
Hall states observed in the lowest Landau level are char- 
acterized by a filling factor with odd denominator may 
be attributed to the fact that a stable cluster made up 
of p electron coordinates and h hole coordinates can only 
be obtained if one of the parameters is odd, thus corre- 
sponding to case (a) or (b). Indeed, for p and q both 
even, the paths obtained by clustering p coordinates Zi 
and h coordinates Wj could be energetically unfavorable 
compared with paths obtained by clustering p/2 coordi- 
nates Zi and h/2 coordinates Wj. 

The incompressible state at v = 5/2 is characterized 
by p = h = 2 ( case (c) ). The,-absence of incompressible 
state at v = 1/2 is attributedtd to the fact that, in the 
lowest Landau level, clusters with p = h = 2 are unstable 
toward the formation of clusters with p = h = 1. Since 
the paths of such p = h = 1 clusters contribute to the 
partition function with a phase (— l)[ p l reminiscent of 
Fermions, they can be regarded as composite Fermionsu. 

It is important to note that the approximation to the 
density matrix provided by function k of Eq. ([lO]) is un- 
able to describe neutral excitations consisting in pairs of 
well-separated fractionally-charged excitations. Indeed, 
when using this approximate form of the density matrix 
to evaluate the number of electrons in a surface large 
compared to the squared magnetic length, one obtains 
fluctuations much too small to include the effect of such 
excitations. In order to include the effect of a single 
quasiparticle-quasihole pair, it is necessary to add to k 
a further term. This additional term can be obtained 
from the right-hand side of Eq. (^) by replacing one of 
the projectors associated to clusters of p electrons and 
h holes by two projectors: One associated to a cluster 
of p' electrons and h! holes and an other associated to 
a cluster of p — p' electrons and h — h! holes, with p' 
and b! satisfying p'h — ph 1 — ±1. The cluster associ- 
ated to the plus sign leads to a local charge defect e/q 
(quasiparticle) and the other to —e/q (quasihole), where 
e is the electron's charge. The parameters p' and h! 
are given by p'/h' = l/(m + l/(n 2 + ... + 1/rik-i))—)) 
in terms of the continued fraction decomposition p/h = 



l/(ni + l/(n 2 + ... + 1/rifc))...)). Thus the creation of a 
neutral excitation corresponds to the breaking of a clus- 
ter of p electrons and h holes. The fact that this break- 
ing requires a finite amount of potential energy, together 
with the fact that k of Eq. (0) considered as a func- 
tion of z\...z p n or of wi...WhN has a small overlap with 
the contribution described above including a neutral ex- 
citation, leads to a gap in the energy spectrum and to 
incompressibility. 

We finally address the question of the temperature 
range in which the function k j (w\...W} 1 nZi...Zpn) of Eq 
.(0) may be adequate to describe the density matrix 
of Eq. (g). Although k is unable to describe neu- 
tral excitations consisting in pairs of well-separated el- 
ementary excitations which are associated to short wave- 
length collective modealil, We do believe that the ef- 
fect of long wavelength collective modes are included in 
K,(wi...WhNZ\...z P N) and are responsible for the correla- 
tions existing between variables Zi and Wj. This collec- 
tive modes are characterized by a finite excitation energy 
£fc=o- We thus conjecture that the temperature for which 
k is adequate to describe the density matrix is of the or- 
der of e k=0 /k B - 
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